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Abstract 



The existence of entangled quantum states gives extra power to quantum computers 
over their classical counterparts. Quantum entanglement shows up qualitatively at the 
level of two qubits. We show that if no entanglement is envolved then whatever one 
can do with qubits can also be done with classical optical systems. We demonstrate 
■".^ ' that the one- and the two-bit Deutsch-Jozsa algorithm does not require entanglement 

' and can be mapped onto a classical optical scheme. It is only for three and more input 

bits that the DJ algorithm requires the implementation of entangling transformations 
and in these cases it is impossible to implement this algorithm classically. 

1 Introduction 

The realisation that quantum-mechanical systems have a very large in-built information 
processing ability and can hence be used to implement computational algorithms has aroused 
much interest recently |^, ^, ^, ^ . The basic unit of quantum information is the quantum 
I bit (qubit), which can be visualised as a quantum two-level system. The implementation of 

quantum logic gates and circuits is based on the tenets of reversible logic and the fact that 
the two states of a qubit can be mapped onto logical and 1 @, |^ . Quantum mechanical 
realisation of logical operations can be used to achieve a computing power far beyond that 
of any classical computer jl^, |ll|, |2[ |l^, Q . A few quantum algorithms have been designed 
and experimentally implemented that perform certain computational tasks exponentially 
faster than their classical counterparts. Three such algorithms are, the Deutsch-Jozsa (DJ) 
algorithm ||l5|| , Shor's quantum factoring algorithm p^ , |l7| , and Grover's rapid search 
algorithm |]18| . 

The existence of entangled states within quantum mechanics is one of the most striking 
features of the theory ||l^. These states have the potential to show nontrivial non-classical 
effects |^,|l|,|2|. It is now well recognised that the EPR paradox is based on the existence 



of entangled states as only such states exhibit correlations that are capable of violating 
Bell's inequalities. It is rather interesting to note that entangled states play a crucial role 
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in quantum computation as well, and it is the entanglement between qubits that gives a 
quantum computer its inherent advantage. A few researchers have discussed quantum algo- 
rithms which do not rely on quantum entanglement and yet are still intrinsically quantum 
in nature; however such computations are not based on qubits p^ . 

We discuss in this paper how entanglement prevents the mapping or realisation of a 
quantum computation using classical waves alone. Consider the polarisation states of a 
classical light beam. These states are in one-to-one correspondence with the states of a 
qubit. All possible states can be realised by using one half- wave and two quarter- wave 
plates. One can thus pass from any desired polarisation state to all others in this way, i.e. 
all U{2) transformations can be implemented using these gadgets [Q. Therefore a single 
qubit has a classical analogue. On the other hand, it is not possible to map all the states of 
a two-qubit system onto the polarisation states of two light beams. The entangled states of 
the two qubits have no classical counterpart. Therefore at the level of two qubits itself the 
possibility of mapping a quantum computer onto classical optical fields breaks down. 

To mathematically define a quantum-entangled state, consider a system composed of 
two parts and described by a density matrix p. If the state of this system can be written 
in terms of the states of its constituent systems in one of the following ways then it is not 
entangled 



In the case (|^), the density matrix p comprises of the tensor product density matrices pi and 
p2, describing the constituent subsystems. Each system has a complete quantum description 
by itself and p is said to be strongly separable. On the other hand, in the case ^ the 
density matrix p though not expressible as a tensor product of subsystem density matrices, 
is a positive sum of such tensor products. The positive coefficients piS can be interpreted as 
probabilities and therefore p can be thought of as a classical mixture of strongly separable 
pieces. Such a density matrix is termed weakly separable. However, if the state p is such 
that it cannot be expressed in either of the forms described in Eqns. ([|) and (^), it is said to 
be entangled. For separable states the correlations present can be given a classical meaning 
by interpreting the positive coefficients piS as probabilities, an interpretation which is not 
possible for entangled states. 

The central theme of this paper is to explore the scope of realising quantum computations 
on classical optical systems. In this context we demonstrate the fact that the DJ algorithm 
for up to two input bits is classical, since an explicit realisation of its implementation on 
a classical optical system based on polarisation is possible. Further we show that for more 
than two qubits the classical optical model fails since the algorithm essentially relics on 
quantum entanglement in this case. We stress that it is only at the level of three or more 
qubits that the DJ algorithm is a "truly quantum" algorithm. 



The DJ algorithm was one of the first algorithms that demonstrated the power of quantum 
computers over classical ones It determines whether a Boolean function / is constant 
or balanced. Classically, the algorithm requires many function calls to solve the problem 
without error. The quantum computer solves the problem using only a single function call. 



p = p^ ® p^ 



(1) 
(2) 



P 




2 The DJ algorithm 
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Consider an rt-bit binary string x] a function / can be defined on this n-bit domain 
space to a 1-bit range space, with the restriction that either the output is the same for 
aU inputs (the function is constant) or the output is for half the inputs and 1 for the 
other half (the function is balanced). All the 2" possible input strings are valid inputs for 
the function (/(x) = {0,1}). In quantum computation, these n-bit logical strings are in 
one-to-one correspondence with the eigenstates of n-qubits, and one can hence label the 
logical string x by the eigenstate \x). Classically, for an n-bit domain space, one needs to 
compute the function at least 2"/^ + 1 times in order to determine whether it is constant 
or balanced. The DJ algorithm achieves this on a quantum computer using only a single 
function call [p^ , [2^ . The original algorithm required the implementation of the function 
f{x) encoded through a unitary transformation on an extra qubit, along with with the 



Hadamard transformation lia, 25 



The Hadamard transformation for n-qubits is given by 

H-\x)=Y,{-lf^^^^''^\y) (3) 

3/=0 

where Xj and yj are the jth entries of the n-bit strings x and y and symbolises addition 
modulo 2. The Hadamard transformation plays an important role in many quantum com- 
putational algorithms and when applied to the state |0)n_bit generates a superposition of all 
possible eigenstates of the n qubit system. For a single qubit, the Hadamard transformation 
reduces to 



|0)^;i5(|0) + |l))^ ^ 1 



|i>-:7?(|0)-|i)) ' ^ 



H = H-' = — [ \ { ) (4) 



The n-bit Hadamard transformation is non-entangling in nature and is just a tensor 
product of n one-bit transformations. 

A modified scheme can be designed to solve the n-bit Deutsch problem, using n qubits 



alone |26|. Here, for every function / a unitary transformation is constructed, such that its 



action on the eigenstates of n-qubits is 

l^)n-bit-^(~l)'^^^l^)n-bit (5) 

Consider n qubits, all in the state |0) ; a Hadamard transformation converts this 
state to a linear superposition of all 2" eigenstates with equal amplitudes and no phase dif- 
ferences. The unitary transformation U / (defined in Eqn. ^) acting on this state, introduces 
an /-dependent phase factor in each eigenstate in the superposition. At this juncture, all 
information about / is encoded in the quantum state of the n qubits. A Hadamard trans- 
formation iJ" is once again applied in order to extract the function's constant or balanced 
nature: 

2"-l 2"-l 



a;=0 x=a 



.T = l/ = 
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The final expression for tlie output state in Eqn. ^ has an amplitude 1 for the state |0)jj^_]-,j^ 
for a constant function and an amplitude for a balanced function. The categorisation of 
the function as constant or balanced through a single function call using n qubits, is shown 
pictorially in Figure 1. The number of functions for the n-bit Deutsch problem is ^ C fq/2~^'^ 




Figure 1: The block diagram for the modified DJ algorithm 

(where N — 2"). The experimental implementation of the modified DJ algorithm for n 
bits requires the realisation of the unitary transformation corresponding to each of these 
functions, and the n-bit Hadamard transformation, on a physical system. There have been 
many experimental implementations of the DJ algorithm and its modified version, using 

NMR 

3 Classical optical implementation of the DJ algorithm 

We now proceed towards analysing this algorithm for 1 and 2 bits in detail, with a view 
to exploring the possibility of realising it on a classical optical system. Let us consider the 
following classical system: a monochromatic light beam propagating in a given direction with 
a pure polarisation. The polarisation states of such a beam are in one-to-one correspondence 
with the states of a two- level system and this beam can therefore be visualised as a qubit. 
It is also well known that all unitary transformations on the polarisation states can be 
performed. Consider a birefringent plate with its thickness adjusted to introduce a phase 
difference of 77 between the x and y components of the electric field and whose slow axis 
makes an angle (j) with the x axis. The unitary operator corresponding to this plate is given 
by 

jj, ,\ / COS0 —sincf) \ f e"'/^ \ / cos0 sin0 \ , . 

U[V,(P)- gjj^^ J Q g-i,,/2 J _sin(/) COS0 ) ^ ' 

For the choice rj = n it becomes a half- wave plate and we denote it by i?^, while for 77 = n/2 
it becomes a quarter- wave plate and we denote it by Q^. It has been shown that all SU(2) 
transformations can be realised on the polarisation states by taking two quarter-wave plates 
and one half-wave plate with suitable choices of angles of their slow axes with the x axis. 
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We will henceforth refer to this device capable of implementing SU{2) transformations as 
Q-H-Q and a detailed discussion is found in Q. 

Further, let us map the x polarisation state to logical 1 and the y polarisation state to 
logical 0. With this mapping we can proceed to work with this system as a qubit. Since this 
system comprises essentially of classical elements, we call it a "classical qubit" . 
One-bit case: For the purpose of implementing the one-bit DJ algorithm on a "classical 
qubit" , we need to realise the Hadamard transformation H which corresponds to an anti- 
clockwise rotation of polarisation by 45° in the x — y plane and the four Uj transformations 
corresponding to the four possible functions. These are certain SU{2) transformations and 
therefore are implementable using the Q-H-Q device discussed above. The sequence of optical 
operations is detailed in Figure 2. As an example, consider the third unitary transformation 

c^r"' = (J \ ) (8) 

which when realised on a "classical qubit" would leave the y polarisation unaffected while 
the X polarisation picks up a phase factor of e*'^. This transformation is achievable by a 
single half- wave plate. 
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Figure 2: Optical realisation of 1-bit DJ algorithm. If the resultant polarisation after the 
second 45° rotation is along y then the function is constant and if it is along x the function 
is balanced. Since these two polarisation states are orthogonal it is possible to distinguish 
them with certainty by placing a y polarizer after the second 45° rotation, followed by an 
intensity measurement. 

The outcome of this exercise is to show that there is nothing quantum about this proce- 
dure, and that the concept of a qubit can be mapped onto the polarisation states of a light 
beam. 

Two-bit case: For implementing the two-bit DJ algorithm, we consider two beams of 
monochromatic light representing two "classical qubits" . The Hadamard transformation 
is again a rotation of the polarisation of both the beams in the anti-clockwise direction 
by an angle of 45°. All the eight [7^"*^'* transformations corresponding to the eight functions 
turn out to be factorisable as a direct product of two SU{2) transformations, one acting on 
each qubit. 

= u}<giU] with U},U] e SU{2) (9) 

More specifically we give in Eqn. |l^ the decomposition of four of the eight functions. The 
other four transformations differ from these four by an overall phase factor of tt and therefore 
can be similarly factorised. 
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Figure 3: Optical realisation of the DJ algorithm for two qubits. Similar to the one-bit case, 
the polarisers along y are placed after the final Hadamard transformation to project out 
the 1 00) state. The presence of light signal in both the beams indicates a constant function 
while the absence of signal in either of the beams indicates a balanced function. 
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(10) 



To implement these transformations on two "classical qubits" requires insertion of a half 
wave plate in one or both beams and can be achieved in a straightforward way. Hence the 
implementation of J7|~^'* can be achieved by applying the required SU{2) transformation on 
each qubit separately. The schematic diagram of this implementation is shown in Figure 3. 
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4 Entanglement at the three-qubit level 



We now turn to the case of three input qubits. In this case too the Hadamard transformation 
corresponds to a 45° rotation of the polarisations of each of the three qubits and can be 
implemented easily. The functions Uj-^^^^ arc 72 in number and any physical implementation 
would require a prescription to implement all of them on a system of three qubits. Consider 
for example, a particular balanced function given by the 8x8 diagonal matrix 
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(11) 



This simple looking SU (8) diagonal matrix cannot be written as a direct product of SU{2) 
matrices acting on individual qubits!. This transformation is entangling in nature and 
possesses the capacity to generate entangled states from non-entangled ones. To see this 
clearly, let us consider the action of this matrix on a simple un-entangled initial state 
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(12) 



The entanglement of this final pure three-qubit state can be demonstrated by computing 
the reduced density matrix for the last two qubits, by taking a partial trace over the first 
qubit 



where(p 
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P 
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(13) 



The result (p^"^)^ 7^ p^~^ implies that the reduced density matrix is mixed, proving the 
entangled nature of the overall state of the 3-qubit system. The question now arises as to how 
do we realise such an entangling transformation for our "classical qubit" system? It turns 
out that there is no way we can realise such transformations on "classical qubits" . It is to be 
noted that there are SU{2)®SU{2)®SU{2) worth of states for the system of three "classical 
qubits" under consideration and this does not contain the entangled states. As shown above, 
entangling transformations on the other hand, generate entangled states from non-entangled 
ones and therefore there is no possibility of being able to construct and implement them on 
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this system of "classical qubits". Hence in order to implement the three bit DJ algorithm 
we need a real quantum-mechanical three-quhit system. Such implementations have been 
discussed in the context of liquid state NMR quantum computing [E^ 0^ 



,30, 31^,11 



5 Concluding Remarks 

A complete mapping of a qubit exists in the classical world of polarisation optics. We 
can consider SU (2) worth of pure polarisation states of a monochromatic beam of light 
as a qubit. All unitary operations required to implement logical operations being SU{2) 
transformations, can be implemented through Q-H-Q. We have shown that the one-bit DJ 
algorithm can be readily implemented on this system. More generally, we can say that every 
operation which can be conceived of for a single qubit, can be executed on this classical 
system. 

We have shown that for more than one qubit the "classical qubit" system can work for 
those computations which do not involve entanglement and for such cases there might be an 
advantage over the ordinary binary computers. In particular for the two-bit DJ algorithm, 
since it does not involve quantum entanglement in its implementation for two bits, we are 
able to realise it on our "classical qubits". This algorithm for the two-bit case solves the 
problem with only one function call as opposed to the ordinary classical algorithm requiring 
three function calls. Hence even when this algorithm allows a realisation based on the 
"classical qubits" it outperforms the algorithm based on ordinary binary logic. Therefore, 
using a pair of "classical qubits" has some advantage! 

Finally we illustrate that an algorithm becomes a true quantum algorithm only when it 
involves quantum entanglement at some stage of its implementation, otherwise it is imple- 
mentable on a set of "classical qubits" . The DJ algorithm for three qubits has been shown to 
involve entangled states for its implementation. Therefore, it becomes impossible to realise 
its implementation using "classical qubits" . 

Acknowledgements: I thank my collaborators Kavita Dorai, Anil Kumar, N. Mukunda 
and R. Simon for useful discussions. 
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